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A so far not considered energy loss mechanism in suspended micro- and nanoresonators due 
to noncontact acoustical energy loss is investigated theoretically. The mechanism consists on the 
conversion of the mechanical energy from the vibratory motion of the resonator into acoustic waves 
on large nearby structures, such as the substrate, due to the coupling between the resonator and those 
structures resulting from the Casimir force acting over the separation gaps. Analytical expressions 
for the resulting quality factor Q for cantilever and bridge micro- and nanoresonators in close 
proximity to an underlying substrate are derived and the relevance of the mechanism is investigated, 
demonstrating its importance when nanometric gaps are involved. 

PACS numbers: 62.25.-g, 62.30.-|-d, 85.85. +j 



I. INTRODUCTION 

High quality factor suspended micro- and nanores- 
onators are required in order to make practical several 
potential applications of such mechanical resonators as 
replacements for electronic filters and reference frequency 
resonators as well as ult rase nsitive mass, force, charge, 
spin, and chemical sensorJ^^l Por microresonators (char- 
acterized by having at least two dimensions in the mi- 
crometer range) the Q values usually range from 0(10**) 
up to O(IO^). For nanoresonators (having at least two 
dimensions in the submicrometer range) it has usually 
been the case that Q hardly exceeds 10'* (Ref. I3|, how- 
ever, more recently, nanomechanical beam resonators set 
to vibrate as nanostrings resulted to have Q ^ A x 10^ 
(Ref. 4) , and nanoresonators based on GaN nanowires vi- 
brating in the megahertz range were reported to achieve 
Q — 4x 10^ (Ref. [5]). As new designs and fabrication pro- 
cesses are created in order to overcome the known energy 
loss mechanism^, specially clamping loss and surface de- 
fects, very high Q micro- and nanoresonators can be ex- 
pected to be available for practical applications. How- 
ever, as known energy loss mechanisms are overcome in- 
creasing the quality factor, new mechanisms previously 
ignored can start to set new limits on Q. 

In this work a so far not considered energy loss mech- 
anism is investigated. This investigation is motivated 
by the fact that in most practical applications the res- 
onators are expected to have their motion driven and 
detected electrostatically, that means capacitively, in de- 
signs involving very small gaps extending over large ar- 
eas between the resonator and the electrodes. For in- 
stance, in the current practical designs of RF MEMS fil- 
ters, sub-100 nm gaps are usually required for adequate 
electromechanical coupling and while gaps as small as 
20 nm where already employecP even smaller gaps were 
envisaged as necessary for MEMS filters operation using 
CMOS drive voltage^. Besides, gaps in the nanometer 
range are a natural consequence of the miniaturization 
toward NEMS filters and other devices. 

The energy loss mechanism analyzed in this work re- 



sults from the coupling between the resonator and the 
nearby structures established across vacuum or air gaps 
by an attractive Casimir force. For instance, in micro- 
and nanoelectromechanical resonators the nearby struc- 
tures could correspond to large area electrodes built on 
top of the substrate and located beneath the resonator, 
as is usually the case for beam resonators, or the elec- 
trodes surrounding a disk resonator^. Due to the cou- 
pling across the gap the motion of the resonator results 
in a time varying force on the surface of the nearby struc- 
tures. This force induces the surface to oscillate at the 
same frequency resulting in acoustic emissions that carry 
away a fraction of the resonator mechanical energy. Such 
noncontact acoustical energy loss was considered previ- 
ously in the context of tip-sample interaction in atomic 
force microscopj'^ due to the van der Waals force be- 
tween metals and is generalized here to the interaction 
between the surface of micro- and nanoresonators with 
its surroundings mediated by the Casimir force, calcu- 
lated using the full Lifshitz theory. In the present work 
this mechanism is analyzed in details for the case of sus- 
pended beam resonators, considered to be located on top 
of a substrate. Because in practice the substrate is much 
larger than the micro- and nanoresonators we consider 
in this analysis it is modeled as a semi-space. As simpli- 
fying assumptions both the beam and the substrate are 
assumed to be made from a homogeneous isotropic ma- 
terial, and their motion is considered adiabatic (purely 
elastic) . 

II. THE ENERGY LOSS MECHANISM 

A. The Casimir force 

The Casimir force, which gives rise to the new energy 
loss mechanism, has the same physical origin as the van 
der Waals force, resulting from the quantum fluctuations 
of the vacuum electromagnetic field^^. However, while 
in a simplified picture the van der Waals force can be 
understood as resulting from the propagation of virtual 
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nonretarded electromagnetic waves, resulting in a short 
range effect, the Casimir force originates from the re- 
tarded waves that act at larger distances, extending the 
range of action of the quantum fluctuations. Because 
the Casimir force becomes relevant in the submicrometer 
range, its impact on the operation of MEMS and NEMS 
has been receiving increasing attentiorfSnH. In general, 
this peculiar force depends on the geometry and the op- 
tical properties of the boundaries, however, our analy- 
sis requires solely the knowledge of the negative (attrac- 
tive) pressure between two semi-spaces as first derived by 
LifshitM The final expression for the force is a function 
of the optical properties of the semi-spaces through the 
frequency dependent complex dielectric function. Using 
the Lifshitz theory the Casimir force between semi-spaces 
made from materials relevant for the fabrication of micro- 
and nanoresonators was calculated in Ref. [M] Follow- 
ing this last work and the references therein, we express 
the Casimir force for real boundaries in terms of a cor- 
rection factor to the pressure predicted for two perfectly 
conducting plates P°(d) = — 7r^?ic/(240(i^), namely 

Pid) = -V{d)^^ = -V{d}^ - vid)P"{d), (1) 

where d denotes the gap between the surfaces, h the 
Planck constant over 2tt, c is the speed of light, and the 
constant C incorporates the constant factors in the above 
expression for later convenience. The factor ri{d) is usu- 
ally referred to as the finite conductivity correction fac- 
tor, derived from the actual dielectric properties of the 
surfaces involved using the Lifshitz theory. For all known 
materials ri{d) < 1, therefore, the pressure between two 
parallel surfaces made from actual materials is always 
smaller than the pressure between perfectly conducting 
plates P'^ . The analysis presented in Ref. [T3] indicates 
the relevance of this correction factor, which is as small 
as 0.088 for silicon surfaces separated by a 10 nm gap, 
and can not be simply ignored. 



B. Acoustic emission and the quality factor 

Here we consider the setup were a rectangular can- 
tilever or bridge resonator of length I , width w and height 
h is placed a distance d above the substrate. When the 
resonator is set to vibrate in a given mode with time vary- 
ing vertical displacement u„(a;,t) = Un{x) expiiojnt) the 
gap varies according to d—Un{x, t) resulting, in the small 
displacement approximation, in a time varying Casimir 
force on the substrate. For an infinitesimal rectangular 
element with length dx this force is 

dF{x,t) = ' wdx 

« C^wdx + 4C^w;u„(x)e^""*da;. (2) 
d^ d° 

The first term represents a constant force and can be ig- 
nored. It is the second time varying term proportional to 
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FIG. 1: Mode shapes for the three lowest frequencies n = 1 
(continuous), 2 (dashed), and 3 (dotted) for (a) cantilever and 
(b) bridge resonators. 



Un{x) which induces a time varying displacement of sub- 
strate surface that, in its turn, results on the emission of 
acoustic waves with frequency aj„. The wavelength A of 
the waves produced on the substrate can be shown to be 
related to the dimensions of the resonator by the approx- 
imate relation A ~ {l/h)l, valid for the frequencies gener- 
ated by the first three modes. Because for most practical 
devices l/h > 10 and w < I, A is large compared with 
the lateral dimensions of the source, therefore, justifying 
the use of the point source approximation. In this ap- 
proximation the details on the force distribution over the 
source are not of fundamental importance for the calcu- 
lation of the irradiated acoustic power. However, some 
aspects of the force distribution must be taken into ac- 
count as we do next. 

Firstly, we note that in most practical applications of 
beam micro- and nanoresonators the ratio h/l is suf- 
ficiently small to the Euler-BernouUi beam theory to 
appljff^, at least approximately. In this case the result- 
ing mode shapes are given by a general expression of the 
form 

Un{x) — UQ{cOsh{KnX/l) — COs(k„x/Z) 

+ Xn [sinh(K„a;/0 - sin(K„a;/Z)]} (3) 

where for the first three modes of the cantilever (bridge) 
we have ki = 1.8751(4.7300), K2 = 4.6941(7.8532), Kg = 
7.8548(10.996), and xi = -0.7341(-0.9825), X2 = 
-1.0185(-1.0008), and X3 = -0.9992(-0.9999). In Fig. 
[Tj Un{x) for the first three modes of cantilever and bridge 
are presented. As seen from the second term in Eq. ^ 
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the amphtude of the time varying force on the surface 
varies in position following Un{x). Therefore, all the 
points of the source region on the surface are in phase 
for the first mode of both the cantilever and bridge, and 
both sources can be held as acoustic monopoles. For the 
second mode of the bridge, points to the left and to the 
right of the mid point vibrate out of phase by 180° and 
the source can be treated as an acoustic dipole. Before 
discussing how the acoustic emissions produced by the 
other vibrating modes can be held, let us consider the 
acoustic emission by the monopoles and dipoles on the 
surface. 

Miller and PursejEH were the first to derive from the 
elasticity theory an expression for the acoustic energy 
emitted by a point source vibrating normal to the surface 
of a semi-space. For an harmonically varying force of the 
form F{t) — F expUcOnt) they employed the admittance 
method obtainingi^ 



1 Vp^i 
47r Cf 2 



F^ 4>M 



where 4>m{l) corresponds to 



0m (7) = -f"! 



■Fo(P,7) 



dp 



with 



(4) 



(5) 



(6) 



and 7 = \/Cli/Cl^, is the density, Cli and C44 are the 
elastic stiffness coefficients that characterize the isotropic 
material of the substrate, as indicated by the superscript 
s. Several authors have rederived Eq. Q since the now 
classical work of Miller and Pursey, h owev er, differying 
in the definitiorfi^ or on the evaluatioriSMI gf 0,„(7). In 
the original worlP^ the integral was evaluated numeri- 
cally taking into account the branch-points p = 1,7, the 
principal value of the radicals, and the only physically 
relevant pole satisfying the condition p > 7. Following 
the prescriptions given by Miller and Pursey in the Sec- 
tion 7 of Ref. I16f a) we derive another representation 
suitable for numerical evaluation 



0m (7) = 



(2p2_^2)2+4p2yfT 

4p3(p2 - 1)^72 -p2 
(2p2 _ ^2)4 16p4(p2 _ 1)(^2 _ p2) 



-dp 



dp 



(7) 



where the last term corresponds to the contribution from 
a clockwise indentation around the pole at p — pr, deter- 
mined as the root of Fo(p, 7) = satisfying pr > 'f. The 
above representation was checked to reproduce the nu- 
merical result reported by Miller and Pursey for 7 = 



4>m{V^) = 0.537, and to differ only sligthly from the re- 
sult reported by HunteJ^^ for 7 = 2, (j)m{2) = 0.415, in 
which case we obtain 0m(2) — 0.409, a difference that 
may be due to numerical precision. It is worth to note 
that in Ref. 17, in spite of the fact that the authors base 
their analysis on the work of Miller and Pursey, the ex- 
pression for the acoustic power 11"* contains an integral 
that, while similar to that in Eq. ([s]), leads to significant 
discrepancies in the numerical results, underestimating 
the emitted acoustic power by as much as a factor of 
10. In a recent analysis of the support (or clamping) 
loss in micromechanical resonatorJi^, which also follows 
the work of Miller and Pursey, exactly the same expres- 
sion for n™ given in Eq. was reported. However, 
the numerical evaluation of 0m (7) did not take into ac- 
count the contribution from the pole at pr, also lead- 
ing to a significant underestimate of acoustic emissions. 
The neglecting of the contribution from the pole corre- 
sponds, physically, to neglecting the contribution of the 
Rayleigh surface waves which are responsible for carry- 
ing away the major fraction of the acoustic energjff^. In 
Ref. ini the acoustic emission by an harmonically vary- 
ing normal force applied to the surface of a semi-space 
was reconsidered. This work provides a coefficient of 
friction Fj^ = (_j_F^ /(Airp^'cf), where q = \fClJ~ff is 
the transverse sound velocity, and f corresponds to the 
sum of three terms similar to those in Eq. ([t]). From 
F_L we can derive the emitted acoustic power defined a^, 
n — Y \2uj^u^ were, in the notation of Ref. f9l wo de- 
notes half the amplitude of the vertical motion which 
was given as the sum of a complex amplitude propor- 
tional to mq plus its conjugate. The resulting expression 
for n™ is identical to that in Eq. Q with 0m (7) re- 
placed by ^4^/(27). By means of an adequate change 
of variables, 0m (7) can be made equal to ^_l/(27) , ex- 
cept for the limits of integration of the two integrals in 
Eq. ([7]). Therefore, the emitted acoustic power calcu- 
lated using the expressions provided in Ref. IHl does not 
match n™ determined by Miller and Pursey. However, 
the results can be made to coincide if the limits defined 
in the two integrals in are taken squared, in which 
case they become the same as the limits in Eq. ([t]) after 
the proper change of variables, indicating that should 
be corrected in this manner. 

Specializing to the case of a suspended resonator vi- 
brating transversally in the mode the time varying 
contribution from the total applied force on the surface 
is 

F{t) = i^exp(iw„t) 

— 4:Cr]{d)d~'^w / Un{x)dx exp{iujnt) , (8) 
Jo 

Therefore, the energy lost per cycle is A[/„ — n//„ = 
2'!Tll/ujn- The resulting quality factor Q is a measure of 
the ratio of the vibrational energy of the resonator C/„ 
to the energy lost, namely, Q — 27rt/„/AC/„ = aj„J7„/n. 
The vibrational energy for both cantilevers and bridges 
is Un = hwlp'^uj'^uo/2 while the mode frequency is = 
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K^E^ h'^ / {12p^ l"^) , implying that for the modes considered 
as acoustic monopoles 



0™ = 



1 



a 



44 



1/2 



wP rj{dY ' 
(9) 

where lu^ = jf^Un{x)dx/(uQl), E"^ denotes the Young 
modulus, and the density of the resonator as indi- 
cated by the superscript r. The result expressed in Eq. 
([9]) reveals that Q"* is a fast varying function of the 
gap distance showing an explicit dependence that goes 
as d}^. However, in order to determine the actual depen- 
dence of Q"* on d we have to take into account the term 
77(d). From the results presented in Ref. [T?] it is gen- 
erally the case for conductors and semiconductors that 
77(d) oc d" with a increasing almost linearly from ap- 
proximately 0.65 for d equal to 15 nm to values close 
to one at 1 nm. Therefore, in this particular range of 
distances, Q™ has an exponent for d varying from 8 up 
to a maximum of approximately 8.7. The dependence 
on the geometrical parameters /i, w, and I indicates that 
the new energy loss mechanism is more relevant for thin, 
wide, and long structures. It can also be inferred that 
is smaller for soft materials in the substrate due to 
the dependence on Cl^ . 

We turn now to the analysis of a dipolar excitation at 
the surface of the substrate. In this case there is no net 
vertical force, instead there is a net bending moment M 
which, due to the symmetry of the second mode for the 
bridge, results to be 



M ^ 8C#u; 
d5 



1/2 



U2{x) 



dx. 



(10) 



This bending moment causes the surface to twist. For 
small sources, we can use the average twisting angle cal- 
culated by BycroflP^ for an harmonically varying bending 
moment M{t) = M exp{iijjnt) distributed over a circular 
region, Eq. (191) of Ref. [501 In the limit of small radius 
over wavelength ratio the amplitude of the angle is 



47r7C| 



where 



dp, 



(11) 



(12) 



and k = 27r/A* = w„ \/ p^ IC\^. We implicitly incorporate 
into the integral the explicit contribution of the physi- 
cally allowed pole at introduced by Bycroft. This last 
author also introduces the same explicit contribution into 
the expression for the average vertical displacement due 
to an harmonically varying normal force when compared 
to the result obtained by Miller and Pursejff^. With this 
definition for the integral we can proceed to obtain a 
representation suitable for numerical evaluation follow- 
ing the same procedure adopted for 0^(7). As we argue 



next, in order to determine the emitted acoustic power 
only the evaluation of the imaginary part of (\)'^(^^) is re- 
quired, which can be written as 



Ap^lp^ - 1)^/72 -p2 



dp 



(2p2 - -y2)4 lQpi(p2 

pI\JpI - 1 



l)(72^p2) 



dp 



(13) 



In the complex notation the real and imaginary com- 
ponents of the displacement correspond to the in-phase 
and out-of-phase components relative to the applied force 
or moment. Only the out-of-phase component of Q con- 
tributes to the time averaged acoustic power through 
n = {Re\M(t)\Re\B{t)\) = Re[M x 9]/2, where ( ) de- 
notes the time average, M the real amplitude from M{t), 
and 9 = iUnO the complex amplitude of the angular 
velocity. We note that we could have used the analo- 
gous definition for the average emitted acoustic power 
for the normal point source obtaining the same result as 
in Eq. Q. In this case H = {Re[F {t)]Re[z{t)\) , where 
z{t) — zexp{iLUnt), and z denotes the complex average 
vertical displacement given by Eq. (129) of Ref. [TCT aV 
In the case of the dipolar source the resulting emitted 
power is 
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s3x 1/2 



P_ 



(14) 



As for n™ it is worth to compare this result for 11'* with 
those found in the literature. Compared to the results 
presented in Ref. [171 also based upon the work of By- 
croft, the same difference is found concerning the def- 
inition of the denominator in the integral 0'(d). Our 
result can also be compared with a derivation for the 
acoustic energy loss due to an AFM tip vibrating par- 
allel to the surface of a plane substrat^. This is an 
analogous situation because the harmonic horizontal dis- 
placement of the force can be interpreted as resulting 
into a time-varying harmonic torque about an horizon- 
tal axis perpendicular to the direction of the tip vibra- 
tion. Considering the definition of the horizontal dis- 
placement given in Ref. '9', which has amplitude 2ito, 
a point force displacing horizontally results in a torque 
M{t) = 2uoF cos{ujnt) — Af cos(a;„i). From the given 
expression for the friction coefficient (see the Appendix 
B of Ref. El) 



Btt pel 



(15) 



the emitted power 11 = r||2a;^itQ results to be the same as 
that in Eq. ( 14) with ipdil) replaced by (7/2)^11 , where 
is an expression similar to Eq. ( 13 ). As for the monopole 
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case, (pdif) can be made to coincide with (7/2)^11 after 
an adequate change of variables except for the hmits of 
integration. The two resuhs for 11'^ can be made identical 
if the limits of integration in ^|| are taken squared. 
From Eq. ( 14 ) we can follow the same procedure as 



for the monopole in order to calculate the quality factor 
for the second mode of the bridge which results to be 



, _ 1.222 2 



(-is r3 



1/2 



1 d^o 
wl rj{dY 



(16) 



Compared to the Q™ calculated for the first mode of the 
bridge, Q'^ is larger by roughly a factor {l/h)'^. This fac- 
tor is exactly what would be expected from the dipolar 
nature of the source. The power irradiated by an acoustic 
dipole where the two sources are separated by a distance 
D = 1/2, as is approximately the case here, is propor- 
tional to {D/X")'^ — (Z/2A'')^ times the energy irradiated 
by a single monopole with the same strengttPil. There- 
fore, as A* cx P/h it results that Q'' oc {l/h)'^Q"^, as 
noted above. In fact, this relation between Q"^ and Q"* 
is generally valid and because in most of the bridge and 
cantilever resonators found in the literature the ratio l/h 
is close to or larger than 10, the energy loss tends to be 
larger for resonators vibrating in such a way as to pro- 
duce a net vertical force on the surface of the substrate 
as compared to a net bending moment. 



C. Beyond acoustic monopoles and dipoles 

The expressions for Q™ and Q"^ as derived above are 
strictly valid for vibrational modes resulting in acoustic 
monopoles and dipoles, respectively. However, we can 
expect that Es. ^ and ( 16 ) provide approximate results 
for slightly more complex vibrations of the resonators 
whenever net vertical forces or net bending moments are 
the prevailing disturbances acting on the surface of the 
substrate. This fact allow us to extend the results for 
some higher order vibrational modes. The estimate of 
Q for higher order modes is important because the use 
of such modes in practical devices is becoming an alter- 
native as a means to achieve high frequency operation, 
specially in the UHF rang^^Hl Focusing on the first three 
modes of the cantilever and bridge resonators, we can 
firstly argue that the second mode of the cantilever is 
going to lose energy predominantly as a monopole due 
to the net vertical force produced on the substrate, while 
the portion of the cantilever vibrating out of phase emits 
energy as a dipole at a much smaller rate. 

In order to clarify this argument, we note that the net 
vertical force given by Eq. ^ is the same for every mode, 
the difference coming from the integral over Un{x). As 
noted after Eq. (|9| this integral can be written as Iu„uqI 
were J„„ = 0.783 for the first mode of the cantilever and 
equal to 0.434 for the second mode. Therefore, the net 
vertical force produced by the second mode is large, com- 
parable to the force for the first mode, resulting in acous- 
tic emissions that exceed any dipolar emissions produced 
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FIG. 2: Contourplot for the quality factor Q as a function 
of the gap d and length I of the resonator for a fixed ratio 
l/w — 5, and thickness h — 0.1 pm. Resonator and substrate 
made from (a) polysilicon and (b) gold; continuous(dashed) 
lines are for bridge (cantilever) resonators; the contours are 
for Q = 10* (bottom curve), Q = 10^ (middle curve), and 
Q — 10^' (upper curve). 



by the small portion of the source close to the free end 
of the cantilever. In its turn, the third mode of the can- 
tilever produces a surface force distribution that is close 
to the dipolar source produced by the second mode of 
the bridge. However, in this case, due to the lack of sym- 
metry of this mode there results both a net vertical force 
and net bending moment. The vertical force is reduced 
compared to the lower order modes being, in this case, 
proportional to — 0.254. The bending moment at 
the more characteristically dipolar portion of the source, 
limited to the left of the second node at a; = 0.868/, is 
close to that found for the second mode of the bridge. As 
a consequence, for two limiting cases where the ratio l/h 
is sufficiently small (large) that the predicted energy loss 
due to the bending of the surface is much larger (smaller) 
than that due to the net vertical displacement the qual- 
ity factor can be estimated using Eq. (161 [Eq. Q]. 
Finally, the third mode of the bridge can be treated ap- 
proximately as an acoustic tripole, a source comprised of 
three monopoles, two inphase and one out of phase by 
180°. In this case the power emitted is that produced by 
a single monopole that causes a net vertical force propor- 
tional to Iu„ = 0.364. 

From the analysis presented so far, some general trends 
on Q for higher order modes can be advanced. In general, 
the acoustic sources at the substrate surface are going to 
be n-poles, corresponding to the n antinodes, each pole 
having approximately the same shape, and consequently 
intensity, along the resonator. For n odd there is a net 
vertical force that decreases significantly for n > 3 re- 
sulting in a very large ratio kn/Iu„ and, therefore, on 
the increase of Q with n [see Eq. For n even, there 
is no net vertical force and an increase of Q with n, pro- 
portional to {D/X^)~^ ~ (l/h)", is expected based on 
general results for the acoustic emission by multipole^^H, 
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III. RESULTS AND CONCLUSIONS 

In order to illustrate the relevance of the new energy 
loss mechanism we present in Fig. [2] contourplots for 
different values of Q. We consider cantilever and bridge 
resonators made from polysilicon and gold. Polysilicon is 
chosen because most MEMS and NEMS are made based 
mostly on this material, and results based on it are rep- 
resentative of other forms of silicon and semiconductors 
like gallium arsenide and germanium, due to their similar 
optical and mechanical properties. Gold is a representa- 
tive of the class of soft (small Young modulus) materials 
that are also employed in MEMS and NEMS. Another 
relevant feature of gold is its high optical reflectivity, 
which results into a stronger Casimir force coinpared to 
a semiconductoi'i^. The results shown in Fig. [5] are for 
resonators with a constant aspect ratio l/w — 5, which 
is large enough to be representative of a wide number of 
practical resonatorJ^i^ but still sufficiently small for the 
point source approximation to apply. Three values of Q 
were chosen, encompassing values that are currently ob- 
tained for micro- and nanoresonators {Q = 10^ and 10^) 
and those expected from technological improvements in 
future devices {Q = 10^). 

What is revealed by Fig. [2] and other similar analysis 
we performed is that the new energy loss mechanism can 
be expected to be more relevant when the gaps involved 
are smaller than approximately 10 nm. Due to the strong 
dependence of Q on the gap distance, this result holds 
also for structures with a thickness considerably larger 
than the one we considered in Fig. [2] (/i = 0.1 /xm), since 
a small decrease in the gap suffices to compensate for 
large changes in h. Because the predicted quality factor 
is smaller for small gaps and thin resonators, this mech- 
anism should be more relevant for nanoresonators actu- 
ated electrostatically, since in this case nanogaps would 
arise naturally. In fact, as mentioned in Sec. |T] electrode- 
to-resonator gaps as small as 20 nm were successfully fab- 
ricated for electrostatic actuation and readout of RF sig- 
nals using a blade nanoresonator. It is worth to mention 
that this nanoresonator was very long (/ = 30 ^m) and 
had a base thickness of about 1 ^m, demonstrating that 
tiny gaps can be produced for comparatively large (order 
o micrometer) structures, and indicating that sub-10 nm 
gaps can be a common feature in near future micro- and 
nanoresonators . 

Due to the fact that the resonator mechanical energy 
is transferred over a vacuum gap and dissipated into a 
nearby structure, this energy loss can be considered as 
an example of noncontact frictiorP^l. Noncontact fric- 
tion has been considered so far mainly in the context 
of noncontact atomic force microscopy (nc-AFMpl, and 
several mechanisms involving the tip-sample interaction 
were considered to be the source of the energy dissipation. 
Acoustic energy loss was considered in this contexiP^l, but 
the results derived so far are valid for small tip oscil- 
lations (harmonic approximation) and can not be com- 
pared with, for instance, the precise experimental data 



presented in Ref. [Ml measured using large amplitude nc- 
AFM. Other possible sources of noncontact friction are 
anelastic processes on the tip and sample, van der Waals 
friction resulting from fluctuating electromagnetic field'^, 
and electrostatic friction involving electromagnetic emis- 
sions or the Joule effect. However, the bulk of experimen- 
tal data can not be consistently explained by any one of 
the known noncontact friction mechanisms^^-^S 26-28^ 
stead specific models (usually purely phenomenological 
models) of energy dissipation were used to explain the 
data for each experiment. However, some of these non- 
contact friction mechanisms could also contribute signif- 
icantly to the total dissipation of micro- and nanores- 
onators. 

An approximate but straightforward comparison be- 
tween some of the different noncontact energy loss mech- 
anisms can be done by comparing the friction coefficient 
per unit of area resulting for each mechanism. This 
coefficient can be obtained by modelling the resonator 
as an one degree of freedom system subject to a vis- 
cous damping. For the first mode of both cantilever 
and bridge resonators the coefficient of friction due to 
the acoustic losses is given by 7 = m^ff uji/{QS), where 
m^f f — c p hwl denotes the effective mass, with c = 0.396 
for bridge and c — 0.250 for cantilever, Q corresponds to 
the quality factor and S = wl to the area. The resulting 
7 is proportional to the area, and for a microresonator 
made from gold with / = 5w — 5/im suspended 10 nm 
above a gold substrate it results to be 7 = 0.32 Kg s~^ 
. For comparison, the corresponding friction coeffi- 
cient due to the van der Waals friction assuming clean 
gold surfaces is approximately 7"'^^ = lO"'^ Kg s~^ m~^ 
at a temperature T — 300 KP^. For semiconductors j'"'^^ 
increases, and for silicon carbide it is predicted to be 
one order of magnitude larger than for good conductors. 
Almost exactly the opposite of that is observed for 7 
which is approximately one order of magnitude smaller 
for a semiconductor like silicon as compared to gold. It 
is worth to mention that the van der Waals friction is 
expected to increase by orders of magnitude under cer- 
tain circumstance^^, for instance, with surface contam- 
ination, therefore giving rise to a significant energy loss 
channel, possibly comparing to or surpassing the energy 
loss mechanism analyzed in this work. It is also inter- 
esting to note that the Joule dissipatiorPsl^ due to the 
Joule effect, investigated in the context of nc-AFM, has 
already been incorporated into the modelling of practi- 
cal micro- and nanoresonators actuated electrostatically. 
In this last case the time varying electric field, due to 
resonator vibration, results into a time varying charge 
at the electrodes and, consequently, an electric current. 
The vibrational energy is dissipated by the Joule effect as 
this current flows through the structure facing the elec- 
tric resistance R forming the equivalent RLC circuiP^. 

We conclude by noting that the energy loss mechanism 
we investigated can be relevant for a wide class of future 
NEMS and MEMS where moving parts are separated by 
distances at the nanoscale. The ubiquitous Casimir force 
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can produce the coupling between the moving parts and 
nearby structures through which mechanical energy can 
be lost in ways that were not addressed in this work. 
Therefore, further investigations on the implications of 
this energy loss mechanism on different systems should 
be performed. 
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